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Digital Computer Laboratory- 
Massachusetts Institute of Technology 
Cambridge 39 » Massachusetts 



SUBJECTS BIWEEKLY REPORT. MAI 2> 1954 

To; Jay W* Forrester 

From: Scientific and Engineering Computation Group 

1. MATHEMATICS, CCDIHG Am APPLICATIONS 

!•! Introduction 

During the period covered by this report 338 coded programs were run 
on the time allocated to the Scientific and Engineering Computation (S&EC) 
Group, These programs represent part of the work that has been carried 
on in kZ of the problems that have been accepted by the S&EG Group* 
Progress on 26 of these problems is given below in terras of programming 
hours, minutes of computer time, and progress reports as submitted by the 
programmers in question* 

This past biweekly period was marked by a significant increase in the 
number of programs that have passed from a testing stage to actual pro- 
duction runs. In particular the revised Comprehensive System (GS II) has 
been operating "successfully. Memoranda describing various aspects of 
the revised system will be made available as soon as possible. One such 
memorandum (M~274l), intended primarily for computer operators, has already 
been distributed. Summaries of the new input program, block recording 
(on magnetic tape) program, and modifications for the post-mortem program 
are included under problem #100« 

Progress reports appear below for the first time for eight new prob- 
lems. Problem #162 (for Dr, F. J, Eppling of the Laboratory for Buclear 
Science) is carrying out a phase shift analysis of the elastic scattering 
of protons by 0-^^, Problem #169 represents a study by E,C,H€y of the DGL 
Staff of the utilization of a general-purpose digital computer in switch- 
ing-circuit design. In problem #182„ a 30th order matrix is being inverted 
for Dr, S, C, Abrahams of the MIT Laboratory for Insultitibn Research* 
Of the other new problems, #1?^ and 175 are being done with the Solid State 
and Molecular Theory Group, #178 with the Instrumentation Laboratory, #179 
with the Aerelastic Laboratory, and #180 with students in the Electrical 
Engineering Department* 
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1«2 Programs and Computer Operation 

The following stomnary is included as a guide for interpreting the 
abbreviations used below. A more detailed description of the terms 
involved can be found in M-2^7« 

a* The upper case letter following the problem number has the 
following significances 

A implies the problem is HOT for academic credit, is Msponsored* 
B implies the problem is for academic credit, is Msponsored. 
G implies the problem is IQ2 for academic credit, IS. sponsored* 
P implies the problem is for academic credit, IS. sponsored. 

The absence of a letter indicates that it is an internal S&EG problem. 

b» DIG denotes the Division of Industrial Cooperation. 
DGL denotes the Digital Computer Laboratory. 
GMMC denotes the Committee on Machine Methods of .Computation. 
DDl denotes the Division of Defense Laboratories. 

100. Comprehensive System of Service Routines , developed by the S&EG 
Group at the Digital Computer Laboratory for the input con- 
version of suitably prepared punched paper tapes. When so 
requested, these routines automatically provide a program with 
suitable programmed arithmetic, cycle-counting, and output 
facilities • 

sDCL Staff s Arden, 40 hours; Best, ^(:i hours; Combelic, 61 hours; 

Demur jian, 17»25 hours; Denman, 3^ hours; Helwig^ 70 hours; 

Kopley, 6 hours; Porter, 2 hours; Siegel, 58 hours; "WWI, 125^ 

minutes 

The S&EC system of utility tapes has been successfully operating 
during the past biweekly period. 

A memorandum, M-274^1, describing the -new Comprehensive System 
(CS II) primarily for computer operators, is now available. 

Similar memoranda directed to tape room personnel and to programmers 
will b e issued shortly. 

Further additions to the system are planned but these will be 
tested independently of the present working system. 

Staff 

A new input program has been recorded on group 11 of the auxiliary 
drum. 

The following kinds of paper tape can be read in automatically on 
the photoelectric tape reader (PETR): 

1. CS II Flexo tapes (identified by an initial fc or a fence) 

2. CS II 'Binary tapes (identified by an initial fb) 
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3* Post-mortem request tapes (identified by an initial fp) 

4» Smnmer session tapes (identified by an initial fs) 

5» 556 tapes (identified by a-ny legal initial character which 
is not an f or a fence) 

The additional programs required for the read in of the various kinds of 
tapes are stored on magnetic tape unit and brought into CM by the input 
program. 

Several other modes of operation can be obtained by inserting 
information into the right insertion register and pressing the upper 
activate button. Among these ares 

1, automatic read in pf paper tapes on the mechanical reader; 
2» recording a stop character on unit 3? 

3» examining a selected register on the magnetic drums; 

4* select a given utility program from unit and read it into 
CM. (These include the scope calibration program, a drum checking program, 
a magnetic tape checking program, and a PETR test program); 

5« obtaining a scope post-mortem; 

6» manual post-mortem requests* 

Automatic logging is provided in the case of paper tape read in« 
For S&EG tapes (whose titles contain dashes) this will take ^lace on the 
scope and the direct punch. For all other tapes (whose titles contain 
m or p) this will take place on the direct typewriter. 

A memorandum describing the detailed operation of the input program 
will be issued shortly* 

Helwig 

A utility program, 100-12-901?, is in use by S&EC staff members 
for doing the following jobss 

1* recording blocks (i.e., parts of the GS U programs, etc.) 
on MT unit in such a form that they can be located and read in by the 
input program. This involves searching for the correct position on 
magnetic tape to begin recording, reading in and optimizing one or more 
55^ paper tapes, and finally recording the equivalent of the paper tapes 
on the magnetic tape as a block with a minimum fof control words. 

2. re-recording blocks on MT unit 0, which differs from recording 
only in that it replaces a block previously a?ecorde<i by a new block which 
can, of course, be no larger than the original block. 

3» optimizing 556 tapes, which consists of reading in one or more 
556 tapes and punching out their equivalent as a single 556 tape which 
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contains no tinnecessary words to be stored. 

This program will read in most tapes that the input program will read in. 
The exceptions are*, (l) the title may contain at most 20 Flexo characters, 
(2) no more than 14 different drxim groups may be addressed, (3) no 55^ 
ditto blocks may occur, (4) at no time during the read in of a tape may 
the words stored form more than 459 groups of consecutive registers* 

The program is being modified to remove restrictions (3) and (4) 
and to permit punching out 556 ditto blocks and to make some minor impirove- 
ments in the way the manual intervention registers are used to control 
the pperation of the program. 

Best 

The following corrections and additions to the OS II post-mortem 
system either have been made or are being tested: 

1. elimination of an automatic stop character after tape post- 
mortem requests; 

2. inclusion of ten carriage returns before any title; 

3. printing a title before the automatic PA PM instead of after 
on fp tape requests (excluding CS I PA FM); 

4« inclusion of an identification mode; 

5. correction of an inactivity alarm which occurred on all 
Manual Intervention Register requests for a GS I PA PM after the first 
fp tape request, 

Arden 

The CS II conversion post-mortem^ which provides information about 
certain kinds of errors detected during the conversion process, has been 
revised, and the latest version recorded on Magnetic Tape unit as part 
of the GS II conversion program* The program apparently operates correctly* 

Siegel, 

An error in the output adaptation section of the CS II conversion 
program was detected and corrected » A great many tests have now been run 
successfully on the output adaptation. 

Demiirjian, Porter 

101 G. Optical Properties of Thin Metal F-h'\m^ on transparent backings 
are determined and printed out automatically by this program; 
the input data consist of the observed reflection and trans- 
mission coefficients 9 the index of the backing, "the waveleng^, 
and the sample thickness. The program calculates by means of 
an iterative procedure and prints out the index of refraction 
and the absorption coefficient of the film, the rate of var- 
iation of these constants with reflection and transmissioi^, 
and the film's conductivity and dielectric constant. 
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sfor Professor L. Harris, Chemistry Department, Dr. A«Ii.Loeb 
5 by Dr. A. L, Loeb(DIG); Eichmond, 45 hours 
iDCLs WWI, 74 minutes 

An "approximate program" had been developed to compute optical 
constants of thin films from reflection and transmission data. This 
program neglects the backing and assumes that the film thickness is 
much less than the wavelength of the radiation. Its results are used 
in the "principal (accurate) program" as first estimates, from which 
the principal program then calculates accurate values. In the visible 
wavelength region the approximation is not very good, but the first 
estimate does not need to be. In the infrared, where the first estimate 
has to be very good, the approximation is excellent. Therefore, this 
combined approximation plus iteration program has proven very satisfactory 
in both the visible and the infrared regions. 

This program is now being "cleaned up", i.e., manual modifications, 
etc. are being incorporated into the program and new printouts are being 
made. Columns and headings are lined up properly, and a provision for 
interrupting the iteration procedure after a preassigned number of iter- 
ations has been perfected. Thus the program and tape will be more amenable 
to general use when the original programmer is not present. 

Large quantities of data obtained in the Chemistry Department on 
gold films are now being run on W»I. Twelve wavelengths in the region 
between 4400 and 6600 i are under consideration. Since the density of 
gold films is not accurately known, four different densities are assumed, 
and the resialts compared. 

106 C . MIT Seismic Pro.1ect is concerned with the development of methods 
for locating deep reflections from underground strata in 
seismic prospecting. The basic method is one of prediction 
by means of an optimum linear operator, 
sfor Professor P.M.Hurley, Geology and Geophysics; Professor 
G. Wadswortii, Mathematics Department 

sby E^A.I^obinson (Res. Assoc.); Briscoe, 36 hoiars; Simpson, 20 
hotirs; Walsh, 25 hours; 
:DCLs Wil, 215 minutes 

The majority of caaputed results during the last biweekly period 
have consisted of spectral and correlation estimates * a continuation of 
a project to categorize the non- stationary properties of seismograms. 
The other two projects mentioned in the last biweekly have suffered delays 
and computations have not progi^^essed. 

In the meantime the group has been training certain of its newer 
members in the use of the programs associated with problem IO6 so they 
will be prepared to carry on experiments during the summer months. 
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113 C» A Stress Analysis of an L«.$faaped Homogeneous Plamar Structure 
is being made for the oase of a concentrated static load. 
The structure is approximated by a framework of bars which 
will deform in the same' Manner as the prototype. This frscae- 
work is then analyzed using the principles of virtual work 
and Southwell relaxation techniques. Boundary conditions have 
been specified for the edge of the framework so that the 
deformations of the model will conform to the actual defor- 
mations of the structure* 

5 for Professor J,SJlrcher, Department of Civil and Sanitary Eng* 

iby S. Sydney (Res, Assist .CMMC) , 20 hours 

sDCLs WI9 121 minutes 

Supplementary programs have been written to obtain data from 
previous calculations in a more usable form for plotting desired graphs • 
Intermediate i^irlwind results and final bar forces are obtained in the 
form of a 55^ tape, facilitating subsequent calculations for stress 
intensities in the prototype plate* Two sets of results have been obtaiB- 
ed with this program* 

119 C, Spherical Wave Propagation produced by the sudden release of a 
spherical distribution of compressed air in the atmosphere 
is being studied by numerical means. This involves replacing 
a set of non-linear hyperbolic partial differential equations 
in 2 independent and 2 dependent variables by a set of difference 
equations written along characteristics. An iterative procedure 
is used to solve these equations* 

:for Professor G* C. Lin,, Mathematics Department 

:by A, Ralston (OfflJC), 5 hours 

sDCLs WI, 31 minutes 

The numerical difficulties mentioned in the past few biweeklies 
have proved more serious than originally thought, Near the shock a 
large percentage error in one very small quantity subsequently causes a 
large error in magnitude in other quantities. Due to these difficulties 
a good deal of reformulation of the problem will have to be done before 
further numerical results can be obtained* 

126 G* A Data Reduction Program for use in the Servonechanisms Laboratory 
is being developed in separate stages to be combined at a 
later date. The first stage is concerned with devising a 
program to fit polynomials to arbitrary empirical functions 
using a least squared error criterion^ The procedure makes 
use of Legendre polynomials and matrix multiplication. 
:for J*E,Ward, Servomechanisms Laboratory, DIG 7I38 AF33( 616)2038 
sby D*T,Ross(DIC) ; Turyn, 65 hours; Hamilton, 10 hours 
sDCL; WWI, 22 minutes 

Part 2 of the main program is operating correctly. The main 
program has been tried once unsuccessfully. Further tests are continuing, 

A scope routine for printing out both single length and (24,6) 
numbers is being written for inclusion in the main program* 
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131 • The Training of Mew Personnel, Tours and Demonstrations are 

among those activities included in this problem • Generally 
speaking, any approved staff problem relating to training 
and/or demonstrations is considered to be in this category^ 
sDCL Staffs Demurjian, 10*5 hours; Kopley, 2 hours; WWi, 72 
minutes 

During this biweekly period computer time was used for demonstration 
purposes and for the filming of a television movie as described below. 

On April 26, twenty members of the Society of Wonen Engineers 
were given computer and Flexowriter demonstrations and a tour of Whirlwind 
installations after a lecture on computers by J. D, Porter • 

On April 27 » nine CBS-TV employees spent the day filming the 
computer for a TV film on MIT research projects to be shown on a program 
called "The Search", A series of such films describing major research 
projects throughout the country will begin* in the fall. 

On April 285 fourteen students from Harvard University were 
given a dfemonstration on and a tour of Whirlwind after a lecture on 
digital computers by C, W, Adams, 

132 C, Subroutines for the Mumerically Controlled Milling Machine are 

being revised and tested. The set of subroutines facilitates 
programming of the computations involved in the preparation 
of numerical data used to control the milling machine. The 
subroutines involve routine numerical and logical operations, 

J for J,0,McDonough5 Servomechanisms Lab,, DIC 6873 

%by J,H,Runyon (E.E, Res Assist), 25 hours 

sDCLs WWI, 12 minutes 

A revised routine for Lagrangian intea?polation which takes 
advantage of extra buffer provisions in CS II was written and success- 
fully tested. Testing of a program -for preparing tape for Series 16 
airfoil templates was begun. Errors which caused the first run to be 
unsuccessful have been found and corrected, 

140, Summer Session System consists of a conversion program, an 

interpretive routine, and mistake diagnostic routines stored 
in Wn^ A special mnemonic instruction code has been developed 
for use with this system, thus simulating a computer with char- 
acteristics quite different from those of WI, This Summer 
Session (SS) computer was developed for the use of students 
participating in the MIT 1953 summer session course on Digital 
Computers and Their Applications, The SS computer was used 
ty the E, E, Department courses 6,535 and 6,25 and is avail- 
able to programmers with suitable problems, 
sDCL Staffs Siegel, 3 hours; ¥WI, 52 minutes 

The SS computer is being used by Mr, R, C. Mills in connection 
with Problem #180, Several SS programs have been run successfully. In 
the course of this work, a mistake in the pat instruction was detected. 
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The mistake was corrected by changing one instruction in the SS PA« 

The SS computer is ready for recording on magnetic tape unit 
as a part of Cs II. It will be recorded as soon as sufficient space on 
the tape unit becomes available* 

1^2 D« A study of Shock Waves has been undertaken in two dimensional 
solids subjected to impulsive loads. The analysis approx- 
imates the solid by two dimensional grid with concentrated 
masses at nodal points* The response of this system is 
computed from a finite difference approximation to the differ- 
ential equations of motion of this system. 

sfor Professor CN^Norris^ Department of Civil and Sanitary Eng. 

sby S, Sydney (Res. Assist, GMMC), 30 hours 

sDCLs WI0 36 minutes 

Shear and stress intensities are being calculated frcm displace- 
ment values of the spring mass system simulating the behavior of an 
extended body under shock loading. The program written for these cal- 
culations has been tested and is being incorporated into the main program, 

Ik? C, Energy Bands in Crystals are being studied by finding solutions 

of the corresponding second order linear differential equation 
satisfying boundary conditions at the origin. The solutions 
are found approximately by using the Gauss-Jackson formula for 
forward integration. The solutions and their first derivatives 
are to be ccanbined in a sum» the wei^ting factors being 
functions of an independent parameter, 

sfor Professor J,G .Slater, Physics Department, DIG 6853 

s by D.J, Howarth (DIC)p 40 hours 

sDCLs WWI, 78 minutes 

Pinal testing of the routine to compute and solve the secular 
equations which represent the last stage of this problem has been com- 
pleted and production work on this routine is now under way* 

155 B, Synoptic Climatology . A multiple regression formula is used 

to predict temperatures from pressure distributions described 
by Tschebycheff polynomials. The matrix of scalar products 
which is used in the calculation of the coefficients of the 
multiple-regression system is being calculated on WWI, 
sfor Professor T.F.Malonej Meteorology Department 
sby R. MillerCDIC)^ 25 hours 
■ sDGL Staffs Porter, 5 hours; ^^0:„ 115 minutes 

Specification of the surface pressure pattern by means of fourteen 
numbers (standardized coefficients of orthogonal polynomials) has been 
completed for the month of January 1948-1952 for a grid the size of the 
eastern half of the United States, 

A program for obtaining the regression coefficients in a least 
squares solution for an arbitrary dependent variable with a given set 
of independent variables has been completed. 
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Crossproducts necessary for determining the predictability of 
pressure 2k hours in advance for 91 points within a grid covering the 
major portion of the North American continent are in the process of being 
computed • 

162 C* Determination of Phase Shifts from Experimental Cross-Sections.. 

A phase shift analysis is being made of the elastic scattering 
of protons by 0-^"« From a partial wave expansion of the 
differential cross-section those phase shifts required to 
provide a fit to a number of experimental angular distributions 
will be determined, 

:for Dr. F.J.Eppling, Laboratory for Huclear Science, Did 6420 

jby E. Campbell(DIC) 9 2k hours 

sDCLs "WWI, 9 minutes 

This problem is concerned with the analysis of a nuclear scattering 
experiment, the elastic scattering of protons by 0^^, over a range of 
bombarding energies from about *5 Mev to 4,6 Mev. Cross-sections were 
measured at eight scattering angles from I68.GO to 90.4®(in the center of 
the mass system) . From the experiment a series of eight curves is obtained 
in which the absolute differential cross-section, d^/dn is plotted as 
a function of the bombarding energy of the incident protons at each of 
the eight scattering angles. From these curves it is then possible to 
obtain d<5'/duj as a function of scattering angle at any given bombarding 
energy* 

For a reaction such as the one above the theoretical expression 
for dd/dus is given as follows s 


















e. 'Si'n 6^ 



JL^I 



where s 



T "T ^ 



= the scattering angle in the center, of mass system, 

Pj( (cos 0) = Legendre polynomial of order j?, 

.J 
P,^ (cos e) ^-^[^—^ Pu(cos e) , and 



M-2812 page 9 

Crossproducts necessary for determining the predictability of 
pressure Zk hours in advance for 91 points within a grid covering the 
major portion of the North American continent are in the process of being 
computed • 

162 C« Determination of Phase Shifts from Experimental Cross-Sections.. 

A phase shift analysis is being made of the elastic scattering 
of protons by 0^", From a partial wave expansion of the 
differential cross-section those phase shifts required to 
provide a fit to a number of experimental angular distributions 
will be determined, 

:for Dr. FoJoEppling, Laboratory for Nuclear Science, Did 6420 

sby E. Campbell(DIC) „ 2k hours 

sDCLs IMLf 9 minutes 

This problem is concerned with the analysis of a nuclear scattering 
experiment, the elastic scattering of protons by 0^^, over a range of 
bombarding energies from about *5 Hev to ^•6 Mev, Cross-sections were 
measured at eight scattering angles from 168,0® to 90«4*^(in the center of 
the mass system) . From the experiment a series of eight curves is obtained 
in which the absolute differential cross-section, dd/dn is plotted as 
a function of the bombarding energy of the incident protons at each of 
the eight scattering angles. From these curves it is then possible to 
obtain d<5"/duj as a function of scattering angle at any given bombardi|)g 
energy. 

For a reaction such as the one above the theoretical expression 
for d5/du> is given as follows s 

ji£> . :X^ -IL esc" ©. e^^^"^"^^ f + 2 CJ^^^^ 5 ^ ^''^ ^^ 

-+-2. i ^ e Str^ d;2 



^»i 



where 



- the scattering angle in the center, of mass system, 
P| (cos 0) = Legendre polynomial of order j?, 
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^ -/ 



|i,r 

j^^ a. ^ - — (t . ^here^j^" is the resonant phase shift and 0^ 

is usually called the potential scattering phase shift or "hard sphere** 
phase shift, ojl is the non-coulomb phase shift of the partial wave of 
orbital angular % and total angular momentum J =J.^l/2; and "where 

' k = 1/x 

'V = the reduced mass of the system, 
V - the velocity of relative motion, 
z,Z = the atomic numbers of the colliding particles, 
£ = the electronic chaiijge 
"h = flank ' s constant divided by 2 IT, and 
/C'^= the orbital angular momentum of the incident protons, 

Using this expression for ^s/d w it is desired to extract tjae o^'^ 
fromil= through J2= 3, i»e., the S-wave, the P-wave, the D-wave, and the 
F-wave phase shifts, explicitly from the experimental angular distributic^ns 
and to determine their {limits compatible with experimental errors, which 
are estimated to be about 15^ to 2^ over most of the energy range covered 
by the experiment. If possible, it is desired to establish the minimum 
number and type of phase shifts required to fit any given angular dis- 
tribution, lie criterion of fit being used is to find the pbikse shifts 
which make the sura of the squares of the percentage errors in the cross- 
section a minimum • 

From these phase shifts it will be possible to determine the 
properties of the excited states of F^? within the region of excitation 
energies which the experiment covers and in addition, to determine, 
within the limits of the experimental errors, whether or not a number of 
liiore subtle factors, such as the ground state of F-^7, a low-lying bound 
state, and possible effects of core excitation, are influencing the 
scattering* 

A preliminary program has been planned to obtain plots of the 
sum of the squares of the percentage errors in the cross-section as a 
function of each phase shift over a range of -ISO© t€3( +I80O, This is 
being done in order to study the behavior of the minimums. So far, one 
tape has been run which gives the error as a function of <4 , the S»wave 
phase shift* 

1^6 G* Construction and Testing of a Delta Wing Flutter Model is being 
effected by replacing the actual wing by a structurally 
equivalent lattice network* An iterative procedure involving 
the evaluation of a matrix equation has been evolved for 
determining the bending and torsional stiffnesses of the 
component members of the network* 

:for M. M. Chen (DIG) 

:by S* Gravitz( Res .Assist. Aero* Eng.), 55 hours 

:DCL Staff: Porter, 6 hours; WWI, I80 minutes 
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Approximately half of the first iteration cycle is now in working 
order. On the basis of the preliminary results » the first trial stiffs. 
nesses were found to yield influence coefficients which matched, on the 
average, to well within 100$^ of the given influence coefficients • The 
second half of the first cycle is completed and is in the process of 
being tested o Testing of this phase of the program has been uneven ' 
because of various machine and programming difficulties. 

The relatively minor modifications to the program to insure 
automatic iteration are being evolved • 

167 Do Transient Effects In Dist411ation are being calculated by 
solving sets of simultaneous non-linear ordinary differ- 
ential equations using numerical integration methods. 
Several of the most iniportant types are being explored 
and results of specific cases will be correla-^ed, 
sfor Professor E« R. Gilliland, Chemical Engineering Department 
sby J.F.O'DonnellCRes.Assist.CMC), 80 hours; Polk, ^ hours; 
Smith, ho hours; Myers, 40 hours 
sDCL Staffs Porter, o5 hour; WWI, 553 minutes 

About half the machine time used by this problem -feras devoted to 
the study of transients in continuous distillations. Some time was 
spent eliminating difficulties in convergence of the steady-state solution* 
¥ith the program and the transient program running successfully, a large 
portion of the data for the thesis of the authors was obtained during 
this period. Also a program was written and operated successfully for 
calculating the change in reflux ratio necessary to maintain constant 
product composition after a step change in feed composition* 

Polk and Smith 

About one- third of the machine time was devoted to the problem 
of take-off in batch distillation. It was spent getting valuable 
exploratory results. This will be continued. 

A program is being written to calculate the same type of results 
as the present program for the limiting case of no-holdup. This should 
be tested in the near future. 

O'Donnell 

The remaining machine time was spent trouble-shooting the program 
for tune-up in batch distillation. Several programming and codiii^' errors 
were eliminated. CS II has been used and some difficulties have developed 
using the output sections. It is hoped that this program will be operated 
successfully soon. 

Jfyers 
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168 D» Indxcial Dowmrash behind a Two-Dimensional Wing , In the analysis 

of the response of an airplane to a sharp-edge gtist, and 
particularly in the calculation of stresses in the horizontal 
tail, it is important to know the downwash at the tail caused 
by the lift response of the wing to the gust. In the present 
solution for the downwash behind a two-dimensional wing, some- 
thing more than the "indicial downwash" is sought. This is 
effected "by allowing the gust front itself to have an arbitrary 
horizontal velocity which, combined with the velocity of the 
airplane, U, results in the wing's penetrating the gust front 
at a velocity Y, 

J for Professor H» Ashley, Aero, Eng, Dept., DIG 672? 

:by N,P,Hobbs (Res, Ass, Aero, lng«), 40 hours 

sDCL Staff? Porter, 1 hour; WWI, 92 minutes 

The quantity CuCA^X*) has been calculated for A+1 X* for four 
values of X*, Physically a this means that the downwash is now known at 
four points in the wake for infinite immersion rate of the wing (? *oo) 
and for times such that the leading edge of the vortex sheet has not 
reached the observation point. It remains to calculate Gjj(A,X*) for 
A+1 X* and to integrate Gjj(A,X*) to obtain the downwash for immersion 
rates other than infinity, 

169 B, Utilizing a General-Purpose Digital Computer in Switching-Gircuit 

Designing to find the most economical physical realization of 
a prescribed switching function. The procedure chosen is one 
that will yield a minimum sum-of -products form . of an^ given 
transmission function within certain limits (i.e,, no more 
than 15 variables or 29 tenas at present), 

sfor Professor D, A. Huff man, Electrical Engineering Department 

sby E,C,Hoy (Res, Assist* E*E,), ^ hours 

sDCLs WWI, 44 minutes 

A minimum sum-of -products form is obtained by applying the 
following rules of Boolean algebra to the functions 

1, .XI + XI« = X 

2. X + X'Y » x+y 
3* X + XI »= X 

4. XI + IZ + ZX' = XI + ZI« 

All but the last of the above rules have been programmed and 
tested successfully. The last has been programmed but is in the process 
of being checked. An output routine has been successfully programmed that 
types out the solution in standard equation form on the delayed printer. 
About ninety percent of the programming has been completed. After the 
last rule is checked out, only a few timing runs need be made to finish 
the problem as far as Whirlwind I is involved. The results will be 
analyzed and incorporated into an SJl. thesis for the Electrical Engineer- 
ing Department, 
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171 C, Improved Power Spectrum Estimates are to be obtained by inves- 
tigating modifications to existing techniques for the numerical 
calculation of Courier cosine transforms to minimize the effects 
of truncation and to supply a confidence curve, based upon the 
interval variations of the calculations, which will help to 
evaluate the significance of the resulting spectrum estimate* 
The object is to obtain a method which will give significantly 
better power spectrum estimates than can be obtained using 
existing techniques* 

sfor J.E.Ward, Project Engineer, Servomechanisms Laboratory 

sby D.T.Hoss(DIC), 80 hours; Hamilton, 20 hours 

sDCLs WI, 163 minutes 

Six functions with known spectra have been run with the Improved 
Fourier Transform program which uses weighting functions^ Most of these 
results will be incorporated directly in an S<,M. thesis report for the 
Electrical Engineering Department <> Two of the runs give the spectral 
windows for the process, which completely describe the properties of the 
methods used. Examination of these results shows that considerable 
information can be gained from the manner in which the successive spectra 
differ. Dense packing of frequencies can be discovered by growth of ■: 
sections of the spectra, for example, A new idea for improving frequency 
resolution through the use of an appropriately chosen "bump" function 
will be briefly tested, although a complete study of this procedure will 
not be possible before the thesis date, 

173. B» Overlap Integrals of Molecular and Crystal Physics » Two-center 
overlap integrals are to be evaluated between various Slater 
atomic orbitals, which are of the form: (power of r) x (exp -ar) 
X (spherical harmonic). By use of prolate spheroidal 
coordinates, formulas in- terms of single functions can be 
derived for these integrals but are unsatisfactory because 
they are of formidable complexity and have false singularities. 
The integrals are evaluated by recasting the expressions in 
terms of spherical Bessel functions of imaginary argument, 
which can be generated by a high-speed computer. 

sfor Professor J. C .Slater, Physics Department 

:by F.J.Corbat(5 (Res, Asst., GMMJ), 80 hours 

sDGLs WWI^ 198 minutes 

The testing of the auxiliary subroutines is now in its final 
stages. In particular, the master-control and one of the sub-control 
routines have been tested. Test results have been obtained for two other 
sub-control routines and have yet to be evaluated. The fourth sub-control 
program is written and typed but untested. Finally the curve fitting 
method was refined and useful fittings of physical data have been achieved. 
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173 • MIT Course 6.537 Spring Term 1954. Twelve students have enrolled 
in the Electrical Engineering course 6»537» entitled Digital 
Computer Applications Practice, being given by Professor 
C« W. Adams, The purpose of the course is to study the advanced 
preparation of coded programs fbr automatic, electronic digital 
computers, in particular, for Whirlwind I, This study will 
include techniques for handling storage and terminal equipment, 
detecting errors and mistakes in programs, the control of 
scale factors, and the use of subroutines* Each student will 
program, prepare on punched tape, and execute on Whirlwind, 
one problem of his own choosing, making use of the CS II 
computer. 
sDCL Staffs Corabelic, 7 hours; Demur jian, 2,5 hours; MWI, 55 
minutes 

#173-100. The problem being programmed is the simultaneous 
solution of a set of linear algebraic equations by the Crout method • 
The upper bound on the order of the coefficient matrix is determined 
only by the capacity of the high speed storage. The solution is carried 
out in floating point arithmetic « The bulk of the program, however, is 
in WWI code, interpretive instructions being used only to form the inner 
products « 

To date, the program has been tested twice on the computer ♦. 
Each test resulted in the discovery of a logical programming error. 
These errors have been corrected and preparations have been made for 
a third trial run on the same fourth-order matrix with known solution. 

Arbuckle 

#173-101. This problem is intended to evaluate the double integral 



'V"=/*ri^<^«] 



cLoJ 

The analysis of the integral requires that the numerical evaluation of 
the error funption be determined. Using an approximation derived by 
the Rand Corporation, tapes have been prepared to determine the precision 
to which this approximation is valid* It is desirable to obtain 8th 
place accuracy. So far, results show values good to 6 places. Initial 
evaluation- of the double integral is being tried for values of /^ = 2 
and a ~ 0. For these values, an analytic solution is readily obtainable 
and is used as a test for the numerical evaluation. 

Brun 

#173-106. This problem is concerned with finding roots of a 
polynomial by Newton's method. The control portion of the routine seems 
to work successfully. However, the roots obtained are not exactly correct. 
At present, the program is in the process of being debugged and more runs 
are being made on the computer. 

Jones 
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#173-107. This problem is concerned with the solution of sets 
of linear equations by the Gauss- Jordan reduction. The program is 
intended to handle dependent as well as independent sets of equations. 
The program has been coded for independent sets of equations, but has 
not as yet given correct results. The routine for solving a dependent 
set has been worked out but has not as yet been coded, 

Kreide 

#173-110 ♦ An Arc Tangent Subroutine has been coded for CS as 
partial credit for course 6,537« By applying Euler's Transformation to 
the Taylor's Series Expansion about zero, the following expression 1?as 
obtained: 



Arc tan x = 



where y = rfe 



ife^ [v,>yAV^*yVv^*..^ 



1+X' 



and V^ = 1 



A V = 2/3 



o 



.n^ = 2 * 4 v6 .., (2|i) 
^ o 3 • 5 • 7 .-. (2n+l) 

This method worked but did not prove as satisfactory as had 
been hoped since about 14 terms are required for x =1 to make the error 
less than 10-5, At present, investigation is being made of the Rand 
Approximations of the Arc tangent.' 

Petrick 

174 C« Tight Binding Calculations in Crystals . The unperturbed wave 

function for a crystal is approximated by a linear combination 
of atomic functions. A Hamiltonian matrix is set up between 
these functions. The eigenvalues of this matrix are the energy 
levels of the electrons. The diagonalization of the Hamiltonian 
will be done on WWI. 

sfor Professor J. C. Slater , Physics Department 

sby G.Koster (DDL), 30 hourS' 

:DCLs ami, 157 minutes 

During the past period it was decided that it was more efficient 
for Whirlwind to calculate the matrix elements of the Hamiltonian matrix. 
These matrix elements are simple polynomials of trigonometric functions. 
This program has been written, tested, and is now being used to calculate 
the energy levels of electrons in nickel. 

175 C. Impurity Levels in Crystals . The calculation of impurity levels 

in crystals involves the solution of difference equations 
which are solved by a Green's function method. This involves 
finding the inverse of matrices (H-E) where H is the unper- 
turbed Hamiltonian and E is the energy of the impurity level. 
The inversion will be done by Whirlwind I. 

sfor Professor J. C .Slater, Physics Department 

sby G. Koster , (DDL) , 30 hours 

:DGLs WWI, 97 minutes 
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A program has been written to calculate the inverse of H-E for 
a series of E values* This program has been tested and found to work* 
It will be applied to the calculation of impurity energy levels in a 
body centered cubic structure. 

178 D« Tra.jectory Study Against an Evading Target . assumes th¥t the target 

evades by flying a circular course instead of the predicted 
linear one. With reference to a certain coordinate system, 
the distances between the missile and the target are computed 
for some initial guess at the independent variable time. The 
rate is observed and when the rate is sufficiently close to 
zero the corresponding value of R is considered the minimum 
passing distance. 

:for Professor H.G.Stever, Aeronautical Engineering Department 

:by C. Block(DIC), 20 hours 

sDCLs WI, 151 minutes 

This is a trajectory study of a missile seeking an evading target. 
What is desired is the minim\am-pas sing-distance of these two objects for 
a variety of missile velocities and target velocities and evasive courses.. 
About 1^0 runs have been made successfully. Another 6OO runs have been 
coaputed by WWI but due to technical programmer's error they are not of 
much use. For each run the minimum- passing-distance, time, 3 parameters 
that define a run, and the number of iterations needed to compute the 
minimum-pas sing-distance are recorded on magnetic tape for delayed printing. 

179 C. Transient Temperature of a Box-Tvpe Beam . The transient temperature 

response of a box beam due to a time varying heat flux input 
on one flange is calculated by an uncoupled finite difference 
procedure. Re-radiation and convection losses are neglected 
and the thermal properties of the structure are assumed constant. 
The results are to be compared with experimental test data. 

sfor J.C.Loria and Professor J.W.Mar, Aeronautical Engineering 
DIG 6727 

sby Lo Ao Schmit Jr^(DIC),,50 hours 

sDGLs WWI, 134 minutes 

The problem is to calculate the transient temperature response in 
a box beam type structure due to a time varying heat flux input on one 
flange. Re-radiation and convection losses are neglected. Variation of 
the local thermal conductivity and specific heat of the structural material 
with temperature is neglected „ Variation of the contact resistance at 
the web flange joint with temperature is neglected. The system does not 
satisfy the thin plate criteria and therefore cannot be treated as a 
unidirectional heat flow problem. The system is treated as two-dimensional 
(i.e., of unit length). The contact resistance of the web flange joint is 
thought to be between certain limits. Transient temperature response 
curves (temperature versus time) have been obtained experimentally for 
five representative points on the structure. The results of these cal- 
culations are intended tos 1.) show which value of the contact resistance 
parameter gives the best correlation with the experimental results? 2.) 
indicate whether the discrepancies arising from the simplifications made 
in this treatment can be tolerated. 
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Tape 179-86-4 j5 ran successfully in 23 minutes (4/3- V54) and gave 
transient temperature response curves at five points in the structxire 
corresponding to the location of thermocouples in the experimental set-up. 




Figure 1 
Schematic Cross Section of Box Beam Showing Location of Thermocouples 

Gomparison of the calculated and experimental transient temperature 
response curves showed that; the calculated results for point A were as 
much as 64<^Fi hij.gher( 1?^ error) than the measured results; the calculated 
results for point B were as much as 36^F lower (21^ error) than the ■ 
measured results; and the calculated results for point C were as much 
as 36^ lower (255^ error) than the measured results. 

At this point more reliable information on the thermal character- 
istics of the structural material became available. This data indicated, 
as did the results discussed above, that a higher value of the thermal 
diffusivity should be used. Adjustment of the specific heat of the 
exposed plate consistant with the new data was also made. Based on these 
changes, which showed promise of an improved correlation, the 179*-86-6 
tape was prepared* 

Tape 179-86-6 ran successfully in 25 minutes (4/23/54) and gave 
transient temperature response ciirves at five points corresponding to 
the location of thermocouples in the experimental set-up» 

Comparison of the calculated and experimental transient temperature 
response curves showed no appreciable improvement in correlation. 

Two factors that help to explain this poor correlation are current- 
ly being studied. The first is that the joint contact admittance between 
the web and flange increases by as much as a factor of three with in- 
creasing mean joint temperature. The detail of the web flange joint and 
the use of steel screws in an aliuninum structure can be shown to produce 
large increases in the joint contact pressure with increasing mean joint 
temperature. 

The second factor previously neglected is that the exposed flange 
should be considered as a thick plate with variable specific heat, rather 
than as a thin plate with constant specific heat. In qualitative terms 
the behavior might be thought of as a protective layer of high specific 
heat being formed thus effectively protecting the inner surface to some 
extent. 

The point at which further use of the canputer will be justified 
depends upon the conclusions reached in re-examining the treatment used 
in the light of the results to date. 
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180 B« Cross correlation of Blast Furnace Input-Outpmt Data. It is 

desired to evaluate a system function relating iron output 
to coke input. This will be done by performing a cross- 
correlation of the quantities and then evaluating the Fourier 
transform of the crosscorrelation function, 

sfor Professor D,P»Campbell and CoA.Powel, Electrical Eng, Dept. 

Jby R. G, Mills, 160 hours; H.J.Scholz,l60 hours 

jDCLs "WWI, 261 minutes 

In an effort to make an analytical study of process efficiency as 
applied to production, the relatively continuous blast furnace process 
was chosen as a vehicle. It is hoped that a transfer function, which will 
relate coke consumed and pig iron produced, Jnay be obtained for the 
process. 

The problem was conceived in an effort to extend the principles 
of servomechanism analysis to the transformation of materials in an 
industrial process. A number of conferences with specialists in blast 
furnace operation has shown that little concrete theory and virtually no 
exact operating techniques have been developed, and consequently, a 
transfer function would be of utmost practical application. Much time 
in the field has been spent to arrive at this result by other methods, 
but no satisfactory conclusion has yet been reached. It is hoped that 
this study will provide a new technique for production process analysis 
and be. of immediate assistance in analyzing the blast furnace operation* 

The method used will be that suggested by Professor X.W.i»ee, 
namely, the evaluation of the Fourier transform of the crosscorrelation 
function relating input and output; these computations will be made 
using 1.) nation-wide pig iron industry data covering a six-year period, 
2,) data for a block of seven furnaces over the same period, and finally 
3.) "much more detailed data for a single furnace over a shorter period 
of time. 

During the past four-week period, preliminary runs, for the piarpose 
of trouble-shooting the several programs, have been completed, and the 
cross-correlation functions for both national and seven-furnace input-output 
data have been computed. Also partially complete is the computation of 
the cross-correlation function for input and output of a single furnace; 
the remaining portion of this computation has been programmed and will be 
run in the immediate future. 

The completed work represents an estimated 75% of the total 
project. Yet to be completed are Fourier transforms of the cross- 
correlation functions. This will be done using a CS program developed 
by D.T.Rgss of .the Servcmechanisms Laboratory and is expected to require 
approximatley 90 minutes of computer time. 

A minor item to be completed is the optimization of a Summer 
Session sine-cosine subroutine. 

Difficulties encountered have been largely programming error, 
although a minor error in the Summer Session computer was discovered* 
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182 C^ Crystal Structures. In solving a crystal structure, the 
corrections to the atomic parameters are obtained by a 
least squares method, as a set of n simultaneous linear 
equations in n unknowns, n may be 30 or more. It is 
required to solve this determinant, 

sfor Dr* .S.C.Abrahams, Laboratory for Ins\ilation Research 

iDCL Staffs Arden, 6 hours; Wfl, 72 minutes 

This problem is a continuation of problem #105 to make use of 
routines developed by Dr. A. Meckler under problem #134, 

th 
The inversion of the 30 order matrix associated with this 

problem is proving difficult because of an eigenvalue spread of about 

1000 ! !• However, it seems probable that at least a crude inverse will 

be obtained ♦ 



1*3 Operating Statistics 

1,31 Computer Time 

The following indicates the distribution of WWI time allocated 
to the S&EC Group, 



Programs 
Conversion 
Magnetic Drum Test 
Magnetic Tape Test 
Scope Calibration 
Demonstrations (#131) 
Total Time Used 
Total Time Assigned 
Usable Time, Percentage 
Humber of Programs 



63 hours, 
20 hours. 



1 hour 



kS minutes 
33 minutes 
40" minutes 
57 minutes 
25 minutes 
12 minutes 



92 hours, 
99 hours, 
93.5^ 
338 



33 minutes 
07 minutes 



1«32 Program Time Distribution 

The following table attempts to show how the WI time expended 
on S&EC programs was distributed with respect to machine runs that gave 
meaningful results (productive computer time) and runs that gave unsatis- 
factory results (lost computer time). Productive computer time is sub- 
divided to indicate the time involved in actual computations as contrasted 
with the time expended getting information out of WWI, Computer time lost 
is subdivided to show the portion of time lost due to errors in the 
programmer's formulation of his problem (logical errors); due to errors 
in the programmer's use of the Wl code, CS Conventions, etc. (technical 
errors); due to tape preparation errors; due to errors by the S&EC computer 
operators in running the program; due to malfunctioning of terminal equip- 
ment; and finally due to miscellaneous causes. 

These times are determined as percentages of the time listed 
above in section 1«.31 ^o^ programs. The times used in computing these 
figures are extracted from the biweekly report forms submitted by the 
various programmers who have used S&EC allocated WWI time. 
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1» Productive Computer Time 

Computation 57*51^ 

Output 13.4$^ 

2. Computer Time Lost Due to Programmers* Errors 

Technical 19.2$^ 

Logical 5.^ 

3» Computer Time Lost Due to Other Difficulties 

Tape Preparation 2.^ 

Operator's Errors 0^ 

Terminal Equipment Malfunction OJ? 

Miscellaneous 1.^ 

1.33 Tape Preparation ( M. Mackey) 

An attempt is being made to obtain some idea of the time expended 
in the preparation of tapes. During the past biweekly period a check was 
made on the tapes processed* 

Due to the variations in procedures involved we have distinguished 
among original complete tapes and the following three types: typed 
modifications - changes of 11 or more registers which must be typed, con- 
verted, then attached to the main program or changes which must be made 
in the body of a Flexowriter tape; manual modifications - changes punched 
directly in 556 form and attached to a converted tape; combined tapes - 
which require duplication of two or more complete tapes. 

The following information was compiled: 



No. of Tapes 

No. of Registers 27901 

Time Constaned 

Thus, it may be seen that the average length of an original complete 
tape is 188.5 registers requiring 34.1 minutes to prepare. A typed 
modification averages 77 >6 registers in length and requires 40.8 minutes 
to prepare while Manual Modifications average 4.3 registers and require 
8.3 minutes for preparation 



Complete 

Tapes 

148 


I^ed 

Mods 

50 


Manual 

Mods 

32 


Combined 

Tapes 

16 


27901 


38?9 


138 




84hrs.8min. 


33hrs.58min. 


2hrs.25min. 


7hrs.8min. 
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2. CCMPUTER EMGBJEERING 

(SoHJDodd) 

Minor troubles have continued to interrupt computer operation 
during the past biweekly period. In many cases the raalfxmctions would 
have been quickly found had they occurred alone, but compounded by wiring 
errors and cold soldered joints, or being of an intermittent nature, they 
ceased 1;rouble incommensurate with their severity. 

A parity check has been installed on the auxiliarj^-storage section 
of the buffer drum. It has revealed an une3cpected number of bad transfers 
from groups ^7» in the form of spurious "ones" JDeing read into the in-out 
register. The source of these spurious pulses has not yet been determined. 

2.1 UMl System Operation 

2.11 Core Memory (A.J.Roberts, L.L.Holmes) 

A considerable amount of time was lost during the past biweekly 
period as a result of sudden tube failures. An intermittent short in a 
buffer amplifier in clock-pulse control which provides synchronizing 
pulses for the in-out system resulted in a^ loss of approximately 12 bourse 
An open filament in a cathode follower in one of the in-out matrices 
caused intermittent program alarms. Several hours were spent locating 
the source of this trouble. As a result of these failures, considerable 
thought is being given to improving in-out trouble- shooting techniques. 
Additional checking facilities may be required to improve the situation. 

Some time was lost because of core-memoiy . parity alarms caused 
by tap shorts in tubes in the digit-plane drivers. 

Loss of power to one of the racjcs was traced to an unsoldered 
connection on an Allied plug-in relay. An investigation of 200 of these 
relays was made. Eight relays were found to have unsoldered connections. 
All the relays of this type in the system will be checked on the next 
installation day. 

2.12 Magnetic Tape (A. X* Perry) 

Reliability of the delayed print-out has been excellent during 
the last biweekly period, much of it due to the marginal-checking facili- 
ties recently installed. Some instances of reported failures were caused 
by programming difficulties. 

All panels necessary for operating Unit 2 as a delayed print- 
out have arrived and will be made ready to operate shortly. Unit 2 has 
been converted to make use of the triangular-hub type reels and is in 
operable condition* S<Me mechanical faults have been observed but are 
not of a serious nature. Efforts to correct these conditions are alrea(^ 
in progress. 
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2*13 Typewriter and Paper Tape (L«H.Norcott) 

Two Flexos were given roiitine periodic inspection and overhaul, 
and the last new Flexowriter was modified for use with V^irlwind during 
the past two weeks. 

Aluminum paper guides have been installed on the back of several 
Flexos* These guides have been affixed with banana plugs so that they 
can be easily removed to permit access to the tape cans, 

2 ♦2 Terminal Equipment (R ,11 ♦Gould) 

Circuits to provide computer control of the output-printer motor 
will be installed soon. Operation will be similar to the control of the 
photoelectric-reader motor. The proper ^ address will turn on the motor 
and insert a delay to let it attain full speed. The motor will remain on 
for a short while after the printer has been deselected and then turn off ♦ ; 
This time will be determined from experience with the autanatic control* 
The operation of the manual switch on the printer will not be affected. 

Modifications to in-out control to simplify the operation of 
block control and the block orders will be postponed for a week because 
of a heavy load of modifications in the rest of the computer. They should 
be done in two weeks. Programs using these orders will not be affected* 

2*21 Magnetic Drums (H*L«2iegler) 

Parity check for groups 4-7, of the buffer drum has been in op- 
eration for the past week and appears to be entirely satisfactory. 

To lessen trouble-shooting difficulties when the fairly large 
electronic he ad- switching system of the auxiliary drym is put into ser- 
vice, a method has been worked out to permit changeover of one digit at 
a time. An additional advantage is that the changeover initially requires 
only one-half the full complement of Type 3 gate-writing amplifiers. 
Because of this modified requirement, changeover of digits began 26 April 
instead of being held up until about 1 June as originally thought necessary* 

An erasing scheme that is simple, sparing of equipment, and easy 
to use has been devised for the magnetic-drum systems. Preliminary tests 
are encouraging. 
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3. AMIl^ISTRATION AMD PERSONNEL 



Staff Transfers 
( J ,0 .Proctor) 

Alexander Vanderburgh has transferred from a DIG staff member 
in 6345 to a DDL staff member* 



New Hon-Staff Personnel 
(R .A. Osborne) 

Jean Harris is a new secretary in the General Engineering Group©* 

Mary Heywood is a new member of the Drafting Department* 

George Kebler is another new member in the Drafting Room* 

Stanley Olsen has joined the Electronic Construction Ship* 

John Pierce is an MIT student working part time in the Tube 
Testing i*aboratory* 

Francis Shaw is a new laboratory assistant in Group 63^5 • 
Terminated Non-Staff 

George Maynard 
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